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ABSTRACT 

Bolometer arrays on large antennas at high, dry sites have unveiled a dusty popula- 
tion of massive, luminous galaxies - submillimetre galaxies, or SMGs - which make 
■ a significant contribution to the star- formation rate density at z > 1. The most cru- 

cial piece of information required to derive the history of obscured star formation is 
the redshift distribution of this galaxy population, N(z), which breaks degeneracies 
in the models and allows the mass and dynamics of the galaxies to be explored via 
high-resolution three-dimensional imaging in CO and by determining their level of 
clustering. Many SMGs are extremely faint, optically; some have no plausible coun- 
' terparts, even in the infrared (IR), making the determination of an unbiased N(z) very 

, difficult. The arrival of the Herschel Space Observatory and next-generation ground- 

based submm cameras will likely exacerbate this so-called 'redshift deadlock'. Here, we 
report the first test of a new method for determining redshifts, based on the observed 
dependence of maser and IR luminosities. We have searched the dusty, lensed, hy- 
perluminous quasar, APM 08279+5255, for the 1612- , 1665- and 1667-MHz hydroxyl 
S— 1 | lines as well as the 22-GHz water line. At z = 3.9 these are shifted to 329, 340 and 

c/2 • 4,538 MHz. Our relatively shallow test data reveal no convincing maser activity but 

we set a meaningful constraint on the OH maser luminosity and we approach the ex- 
pected thermal noise levels, meaning progress is possible. As an aside, we present deep 
new submm and radio imaging of this field. Using a simple shift-and-add technique 
we uncover a new submm galaxy, conceivably at the redshift of APM 08279+5255. 
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1 INTRODUCTION 

The SMG population, discovered using the SCUBA bolome- 
ter array (Holland et al. 1999), has led astronomers a merry 
dance for almost a decade (Blain et al. 2002). Early indi- 
cations, based on the SCUBA Lens Survey (Smail et al. 
2002), painted SMGs as an extraordinarily diverse popu- 
lation, some red, some blue, some starbursts, some with 
unequivocal AGN features (Ivison et al. 1998, 2000). How- 
ever, over the years some common characteristics have be- 
come familiar: on the whole, SMGs are luminous (Lfir ^ 
10 i2-i3 h -2 L ^ wnere hn = # /7i kms^Mpc" 1 ) optically 
faint, red, morphologically disturbed, massive and distant 
(Smail et al. 1999; Gear et al. 2000; Chapman et al. 2003; 
Webb et al. 2003; Neri et al. 2003; Blain et al. 2004; Pope 
et al. 2005; Greve et al. 2005). 

Due in large part to the spectroscopic survey of Chap- 
man et al. (2005) we have seen a good deal of recent progress 
towards determining one of the most crucial parameters in 



galaxy-formation models: the redshift distribution of SMGs, 
N(z). Knowledge of N(z) breaks degeneracies in the models 
(e.g. Blain et al. 1999); it also provides an estimate of their 
typical mass via the clustering characteristics of the popu- 
lation and, for individual galaxies, via observations of CO, 
C I, etc. (Frayer et al. 1998, 1999; Neri et al. 2003; Greve et 
al. 2005; Papadopoulos 2006). 

However, the available spectroscopic redshift distribu- 
tion for SMGs is biased: a radio detection is required to 
provide positional information (e.g. Ivison et al. 2002), and 
rest-frame ultraviolet emission lines are usually required to 
signal redshifts. It is not clear that all SMGs will display 
these traits, certainly not the dustiest systems at the high- 
est redshifts (e.g. Berger et al. 2006). Furthermore, it is not 
clear how redshifts will be determined for the many thou- 
sands of dusty galaxies expected to be detected in wide-field 
surveys planned for SCUBA-2 (Audley et al. 2004) and the 
Herschel Space Observatory (Pilbratt 2004). 

Powerful OH masing is relatively common amongst IR- 
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luminous galaxies: the most recent survey (Darling & Gio- 
vanelli 2002) found that at least a third of ultraluminous 
IR galaxies (ULIRGs, Lfir > 1O 12 L ) support megamasers 
or gigamasers (Loh ^ 10L Q or Loh > 1O 3 L0, respec- 
tively). If starbursts are responsible for a significant frac- 
tion of the luminosity of ultraluminous and hyperluminous 
IR galaxies, as currently thought (e.g. Farrah et al. 2002), 
then the associated turbulence may enable low-gain unsat- 
urated masing (Burdyuzha & Komberg 1990). The earli- 
est OH maser-line observations appeared to demonstrate 
a quadratic relationship between OH and far-IR (FIR) lu- 
minosities, Loh and Lfir. This is believed to be due to 
the abundant flux of FIR photons pumping an OH popula- 
tion inversion in the star-forming molecular gas (Baan 1985, 
1989). Since FIR and radio luminosities (L ra dio) are well cor- 
related (e.g. Helou, Soifer & Rowan- Robinson 1986), emis- 
sion stimulated by the background radio continuum would 
then yield Loh oc LpiRL ra dio oc L FIH . 

Townsend et al. (2001) argued that this quadratic de- 
pendence would yield OH masers detectable amongst the 
SMG population, with peak luminosity densities up to two 
orders of magnitude greater than those seen in current sam- 
ples. This would allow us to determine an accurate and 
relatively unbiased N(z) for SMGs, as well as constrain- 
ing the mass of their black holes, determining their geo- 
metric distances, even probing the evolution of fundamen- 
tal constants (e.g. Barvainis & Antonucci 2005; Lo 2005; 
Kanekar et al. 2005; Caproni, Abraham & Mosquera Cuesta 
2006). Correcting large maser samples for Malmquist bias 
(Kandalian 1996) favours a weaker dependence of Loh on 
Lfir than originally thought, Loh oc Lp'j 2 ^ ' 1 (Darling & 
Giovanelli 2002); on the other hand, using a smaller, com- 
plete sample continues to suggest a quadratic relationship, 
Loh oc Lp' 3 ^ ' 6 (Klockner 2004). Either way, it is reason- 
able to expect OH maser emission from a high proportion of 
SMGs, with Loh in the range 10 4_5 L© for Lfir ~ 1O 13 L 
and mjy-level peak flux densities. 

Maser searches have clear advantages over other meth- 
ods of determining N(z) for SMGs: i) the bandwidth require- 
ment for blind detection of OH megamasers is low, <1 GHz 
(^obs = 165-835 MHz for z = 1-10), smaller still with ad- 
ditional redshift constraints, and the technique avoids con- 
tamination by H I emission from local galaxies (Briggs 1998) 
since there is no significant SMG population at z < 1 (Chap- 
man et al. 2005); ii) the instantaneous survey area is limited 
only by the primary beam (several square degrees for an OH 
line search with the Giant Metre-wave Radio Telescope - 
GMRT) ; iii) interferometry permits some rejection of radio- 
frequency interference (RFI); iv) the position of an emission 
line can be pinpointed accurately, tying an emission line to 
an SMG unequivocably (to a few arcsec, at worst); finally, 
v) the dual-line 1665-/1667-MHz OH spectral signature can 
act as a check on the line identification and the reality of 
detections. 

Here, we report a search for such maser emission. Blind 
searches for maser emission from SMGs will not be feasi- 
ble until next-generation correlators are available, with their 
large instantaneous bandwidths. To begin testing this tech- 
nique we must therefore choose a well-studied, luminous FIR 
galaxy, one with its redshift accurately determined (see also 
Klockner 2004). 

APM 08279+5255 (Irwin et al. 1998) is an obvious can- 
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Figure 1. The IR-to-radio SED of APM 08279+5255, with ra- 
dio data from this paper and a compilation of data shortward of 
5 GHz (Lewis et al. 1998; Dowries et al. 1999; Egami et al. 2000; 
Papadopoulos et al. 2001; Barvainis & Ivison 2002; Soifer et al. 
2004; Beelen et al. 2006). 

didate: a dusty, broad-absorption-line (BAL) quasar embed- 
ded in a gas-rich starburst with a prodigious FIR luminos- 
ity 10 14 Lq). HCN and several CO rotational transitions 
have been detected (Downes et al. 1999; Papadopoulos et 
al. 2001; Wagg et al. 2005) yielding an accurate systemic 
redshift: 3.91126 ± 0.00016 (Axel Weiss, private communi- 
cation). Finally, adding to its appeal, it is amplified grav- 
itationally by a factor A ~ 7 (Lewis et al. 2002) to yield 
an apparent luminosity of ~5 x 1O 15 L0. At z — 3.9, its OH 
and H2O maser lines are both accessible to receivers at the 
National Radio Astronomy Observatory's Very Large Array 
(VLA). Even adopting Loh oc L fi 2 r and an extreme lensing 
model (A — 1000, for which we have no evidence) we would 
expect Loh ~ 3 x 1O 3 L . For a line width of 150 kms -1 
we would then predict a flux density of almost 50 mjy after 
amplification, within reach of current facilities. 

We also report new submm, radio and IR imaging which 
led to the detection and characterisation of a new SMG in 
the vicinity of APM 08279+5255. 



2 SEARCHING FOR OH AND H 2 MASER 
EMISSION 

The choice of configuration for our VLA P-band observa- 
tions (~340MHz) was dictated by the need to reduce con- 
fusion and minimise RFI. We therefore opted for A configu- 
ration, obtaining 4.7 hr of integration during 2002 April on 
the redshifted OH lines. Approximately 80 per cent of the 
data were unaffected by RFI. 

Data were taken with correlator mode 4, with one IF 
pair centred on the 1612.231-MHz hyperfine OH line and 
another IF pair covering the 1665.4018- and 1667.359-MHz 
lines. Each IF had 61 x 48.83-kHz channels, yielding a veloc- 
ity coverage of 2,500 km s -1 . Calibration scans of 0834+555, 
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0542+498 and 3C 274 (1230+123) were obtained, to cali- 
brate phase, amplitude and absolute flux density, and to 
determine bandpass corrections. 

The APM 08279+5255 field is relatively devoid of bright 
radio sources, so we were able to approach the thermal 
noise, with an r.m.s. level of 0.5 mjy beam" 1 in a pseudo- 
continuum image made using all the spectral-line data, and 
as little as 3.4 mjy beam -1 in individual channel maps. In- 
deed, APM 08279+5255 is clearly detected in the pseudo- 
continuum P-band map, with an integrated 334-MHz flux 
density of 5.1 ± 0.8 mjy. It is not clear why, but one IF pair 
produced 2x deeper images than the other. 

We also obtained approximately 10 hr of data with the 
VLA in its C configuration during a flexibly scheduled block 
of time in 2001 August, using the C-band (5-GHz) re- 
ceivers centred on the redshifted 22,235.08-MHz 6ie — 523 
H2O line. Correlator mode 2AD was used: a single IF pair 
with 31 x 390.63-kHz channels giving 760kms _1 of veloc- 
ity coverage. 3C 48 (0137+331) was used to set the flux 
density scale and 0834+555 (5.7 Jy at 4,527 MHz) was once 
again used as the local amplitude/phase calibrator. Again, 
APM 08279+5255 is detected in the pseudo-continuum map 
with an integrated flux density of 1.07 ± 0.10 mjy. Fig. 1 
shows the overall radio spectral energy distribution (SED) 
of APM 08279+5255, adequately characterised by a power 
law of the form S u oc v~ ' 6 , together with submm and IR 
data from the literature. 

For both the C (H2O) and P (OH) spectral-line data, 
image cubes were made using AIPS imagr, with clean 
boxes placed according to the positions of bright sources in 
the respective pseudo-continuum maps. Cubes were made 
with 1- and 2-channel averages and were inspected visu- 
ally for signs of emission. Spectra were extracted as follows: 
APM 08279+5255 and three bright (>30-mJy), compact, 
nearby sources were fitted (jmfit) in the pseudo-continuum 
images with single Gaussians, their sizes fixed to that of the 
synthesised beam (~5.5 arcsec fwhm). Further fits were then 
made in each spectral plane of each image cube, now with 
both object sizes and pseudo-continuum positions fixed. The 
best-fit peak flux densities were thereby extracted for each 
channel. This technique delivers spectra with noise levels 
close to those measured in individual image planes and re- 
mains relatively free of other biases. The resulting spectra 
of APM 08279+5255 were divided by the respective data 
for the three bright, nearby sources, which had been co- 
added and normalised by their total mean flux density. This 
'flat field' was also used to provide a check on the reality of 
any features in the spectra of APM 08279+5255. Noise esti- 
mates in each channel were obtained from source-free regions 
around the quasar using imstat. A simpler approach using 
POSSM gave similar results. 

The resulting spectra are shown in Fig. 2 where no con- 
vincing emission from either OH or H2O can be seen. The 
two weak peaks seen at the expected frequencies of OH (1665 
and 1667 MHz) are no more significant than the dip seen in 
the H2O spectrum at ~100kms _1 . The r.m.s. noise levels in 
the regions around the OH lines at 1612 and 1665/1667 MHz 
are 2.9 and 5.8 mjy, with 0.17mJy for the region around the 
H2O line. These figures compare reasonably well with the 
theoretical expectations: ~2.9 mjy in P band and ~0.16 mjy 
in C band, for good weather and 26 available antennas. 

There is a very weak positive feature at zero velocity in 
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Figure 2. Unsmoothed spectra covering the 1612-, 1665- and 
1667-MHz OH lines (top and middle) and the 22,235-MHz H 2 
line (lower). Uncertainties (l<r) are illustrated with error bars. 
Line positions are marked and velocities are given for an assumed 
redshift of 3.91126. Baselines have not been subtracted. The result 
of extracting a spectrum with POSSM are shown in red in the 
middle panel. 
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Figure 3. Contour map of log 10 Sv, the peak observed flux den- 
sity of an OH maser, as a function of redshift, z, and peak rest- 
frame luminosity density, L v i , adapted from Townscnd ct al. 
(2001), for Q m = 0.27, Q A = 0.73, H = 71 kms -1 Mpc -1 . H 2 
and OH masers, in AGN and IR-luminous galaxies, respectively, 
are shown as triangles and circles, respectively. 



the amplitude spectrum of the cross-correlation between the 
1665-/ 1667-MHz spectrum shown in Fig. 2 and a matched 
filter corresponding to two Gaussian profiles with a separa- 
tion fixed to that of the 1665- and 1667-MHz OH lines. We 
regard this as at best tentative evidence of emission. 

Given the very limited total bandwidth (and the poor 
resulting velocity coverage) we cannot be sure that a 
broad line does not cover the entire spectral window: the 
IRAS 14070+0525 OH gigamaser (z = 0.265, Baan et al. 
1992) would certainly do so, though we would have had 
no problem with a more typical linewidth (~150kms _1 , 
Darling & Giovanelli 2002). However, the peak flux den- 
sity of such an OH line cannot exceed our measured con- 
tinuum levels, which are well below the limits we would 
quote based on the OH spectral noise. For the H2O line, 
where this is not true, we add the continuum level to the 
limit. The peak line flux densities from APM 08279+5255 
are thus 3a < 8.7 mjy for the satellite line of OH at 
1612 MHz, 17.4 mjy for the main lines at 1665/1667 MHz 
and 1.6 mjy for H2O. For a rectangular line of width 
150 km s~ we have constrained the lens- amplified OH lu- 
minosity to below 1.3 x 10 L0 (3<j). As outlined in §1, the 
apparent FIR luminosity of APM 08279+5255 exceeds that 
of IRAS 14070+0525 by three orders of magnitude so a de- 
tection of APM 08279+5255 was quite conceivable and the 
limit we have set for the main lines is significant, lying below 
our most pessimistic prediction. We note, however, that at 
least half of local ULIRGs are not detected and that FIR 
luminosity is clearly not the sole criterion for OH masing. 

Following Townsend et al., the observed flux density S v 
at a frequency v = v 1 / (1 + z) of a maser emitting with a peak 
rest-frame isotropic luminosity density of L v i is given by 
S v = (1 + z)L v i /{^■kD 2 l (z)), where Dl(z) is the luminosity 



distance. The (1+z) factor in this expression accounts for the 
line width narrowing due to the redshift, and partially offsets 
the quadratic drop-off in S v with distance. As an example, 
the peak flux density of lOmJy for the IRAS 14070+0525 
OH gigamaser at z = 0.265 (D L (z) = 4.1 x 10 25 m for 
Q. m = 0.27, Q.a = 0.73, H = 71 kms -1 Mpc -1 ), translates 
into L v , = 1.7 x 10 30 WMHz~\ 

Our limits on maser-line peak flux densities are there- 
fore L v , < 2.7 x 10 32 , 5.4 x 10 32 and 5.0 x 10 31 WMHz" 1 
for the OH lines at 1612 and 1665/1667 MHz and the H 2 
line, respectively. These are shown in Fig. 3 alongside obser- 
vations of masers in other IR-luminous galaxies. 



3 ANCILLIARY DATA: SUBMM, RADIO AND 
IR 

The APM 08279+5255 field was imaged with the SCUBA 
bolometer array (Holland et al. 1999) during 2000 March to 
2001 March, primarily to help characterise the point spread 
function for the deep submm survey of high-redshift radio 
galaxies (M00AU09) presented by Stevens et al. (2003). 

The James Clerk Maxwell Telescope's (JCMT's) sec- 
ondary mirror was chopped and nodded by 30 arcsec in R.A., 
resulting in the characteristic — 0.5/+1.0/— 0.5 beam pro- 
file visible in Fig. 4. After overheads, 25.6 ks were spent in- 
tegrating on the field, split into approximately hour-long 
chunks separated by checks on focus, pointing accuracy and 
atmospheric opacity. The data were reduced using stan- 
dard software, but were imaged using several new tech- 
niques. First, the high flux density of APM 08279+5255 al- 
lowed us to accurately align each 2,560-s of data, peaking 
up at 450 /im having adopted the 1,400-MHz radio position 
for APM 08279+5255 (see later; R.A. 08 h 31 m 41 s .708, Dec. 
+52°45'17.44" J2000). Despite the small shifts being ap- 
plied (1.6 and 0.8 arcsec r.m.s. in R.A. and Dec, consistent 
with normal pointing drifts), this 'shift-and-add' approach 
made an appreciable difference to the final image. Second, 
we employed new versions of the SURF setbolwt and rebin 
tasks, developed for the radio galaxy survey. These enabled 
us to determine accurately weighted and calibrated signal 
and noise values for square pixels of arbitrary size, indepen- 
dently of neighbouring pixels, at 13-arcsec fwhm resolution 
for 850 fj,m. The signal and noise for each pixel correspond 
to the stream of data collected only when bolometers arc 
centred within that pixel in a manner reminiscent of the 
'zero-footprint' approach used by Scott et al. (2002). For 
Fig. 4 we smoothed with 6- and 3-arcsec fwhm Gaussians 
and the resulting 850- and 450-/im maps have noise levels of 
1.7 and 7mJybeam~ (before cleaning). 

This led to the detection of an new, relatively bright 
SMG in the vicinity of APM 08279+5255, with R.A. 
08 h 31 m 35 s .46 (35 s .48), Dec. +52°44'22.4" (26.8") J2000 at 
850 /im (450 /im), with a positional uncertainty of 2.2 arcsec 
at 850 /im (lu, Ivison et al. 2005). We denote this source 
as SMM J083135. 5+524425, adopting the mean submm po- 
sition and following the IAU convention. Using the quasar 
to align the 450- and 850-/im coordinate frames would re- 
sult in a marginal (0.9-arcsec) reduction of the 4.4-arcsec 
offset between the 450- and 850-/im positions. The peak 850- 
/iin flux density of the new source, bootstrapping the cali- 
bration from Uranus (via APM 08279+5255: 286 ± 10 and 
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Figure 4. Grcyscalcs of the APM 08279+5255 field, as im- 
aged by SCUBA at 850 (im (top) and 450 fim (middle) and by 
the VLA at 1,400 MHz (lower). In the submm, contours are 
plotted at —3,3,4. ..10,20. ..100 xa. These are signal-to-noise im- 
ages, smoothed to 15 and 7.5arcsec FWHM. The new SMG, 
SMMJ083135. 5+524425, can be seen to the south west of 
APM 08279+5255, the brightest object. Their -0.5/+1.0/-0.5 
beam profiles are characteristic of real sources in unclcancd jig- 
gle maps. In the radio, contours are plotted at —4, 4... 10 Xcr 
where a is 28 /ijy beam -1 and the synthesised beam measures 
1.45 X 1.16arcscc FWHM with its major axis at position an- 
gle 52°. Circles (8-arcsec radius) mark the 850- fim positions of 
APM 08279+5255 and the new SMG. 



82.7 + 2.2 mjy at 450 and 850 /xm; cf. 84 and 285 mjy, Bar- 
vainis & Ivison 2002), is 11.2+1.7 mjy, making it a relatively 
bright SMG (we would expect around one per lOOarcmin 2 ). 
The flux density rises only marginally, to 11.4 mjy, when 
measured in an 14-arcsec-radius aperture, so any evidence 
that the source is extended is marginal at best. The SMG 
is also detected at 450 fixn with a flux density of 35 ± 7 mjy, 
though only following our shift-and-add treatment. The 450- 
/850-Atm flux ratios of the new SMG and APM 08279+5255 
are consistent (~8) to within the large photometric uncer- 
tainties. 

APM 08279+5255 was then targeted at 1,400 MHz in 
an attempt to further constrain the position and redshift of 
the new SMG. For this we used the VLA in its A configu- 
ration, with correlator mode 4 (28 x 3.125-MHz channels, 
half each for left and right circular polarisations), during 
2001 January. A noise level of 28/xJybeam -1 was achieved. 
SMM J083135. 5+524425 has a single, bright and relatively 
compact radio counterpart centred at R.A. 08 h 31 m 35 s .59, 
Dec. +52°44'27.8" J2000 (~0.2-arcscc uncertainity, la), 
with a total flux density of 386 ± 55 /iJy (after correcting 
for bandwidth smearing and the primary beam response). 

The total 1,400-MHz flux density of APM 08279+5255 
was measured to be 3.05 ± 0.07 mjy so, as with the submm 
flux ratios discussed earlier, the submm/radio flux ratios 
(~28) of the new SMG and APM 08279+5255 are consistent 
with one another to within the photometric uncertainties. 
This suggests they may well lie at a common redshift. A full 
track with the IRAM Plateau de Bure Interferometer, tuned 
to CO(4-3) at z = 3.911, would confirm or refute this possi- 
bility. Although perfectly consistent with a point source, we 
note that the morphology of the new SMG is suggestive of 
an arc and so we cannot rule out the possibility that this is 
a new lens feature caused by an unknown, unseen, massive 
foreground cluster, though we would then expect to see arcs 
from the more numerous optical/IR population. 

We have inspected a Jf s -band image, obtained by us 
during 2000 November 12 in 1.0-arcsec seeing using the 
INGRID camera on the 4.2-m William Herschel Telescope 
(WHT), at the position of the radio emission. The region 
is blank in the near-IR, with K s > 20.5 (3a, 1.5-arcsec 
radius aperture, Vega). Fig. 4 also reveals faint (3<r), co- 
incident 450- and 850-/xm emission, ~15arcsec north east 
of SMM J083135. 5+524425, this time with no radio or IR 
counterpart. 



4 CONCLUDING REMARKS 

We have conducted a preliminary search for OH and H2O 
gigamaser emission from the distant, luminous quasar, 
APM 08279+5255. Our most pessimistic prediction - an OH 
maser exhibiting behaviour of the form Loh oc Lfiri sud_ 
ject to extreme lens amplification - would have been within 
range of our search. We have thus set a meaningful con- 
straint on the OH maser luminosity of APM 08279+5255, 
though we acknowledge that FIR luminosity is not the only 
criterion for OH masing. 

We note that 10 x more bandwidth is already available 
at the GMRT than was used here. Upcoming facilities will 
offer further improvements and we suggest that the Square 
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Kilometre Array 1 , for example, should be designed to have 
sufficient instantaneous bandwidth to employ this technique 
effectively on blank-field SMGs of the type expected to be 
found in their thousands in the coming few years. 

Our submm imaging has uncovered a new, bright SMG 
~75 arcsec from the quasar, possibly within a physical dis- 
tance of <lMpc since their submm/radio flux ratios are a 
good match. Its position has been determined accurately in 
the radio, but it has no IR counterpart. 

The effectiveness of the shift-and-add technique sug- 
gests deep archival searches could be made in the vicinity 
of other bright sources, e.g. JCMT pointing sources such as 
3C 345 and BL Lac that have accumulated many hours of 
integration. Blank-field surveys with the next generation of 
sensitive submm cameras may also benefit from variants of 
this technique, particularly at 450 fj,m. 
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